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SUMMARY 

In this study an attempt is made to put into perspective the problem of a rotating disk, be it a single disk or a 
number of concentric disks forming a unit. An analytical model capable of performing an elastic stress analysis 
for single/multiple, annular/solid, anisotropic/isotropic disk systems, subjected to both pressure surface tractions, 
body forces (in the form of temperature-changes and rotation fields) and interfacial misfits is derived and discussed. 
Results of an extensive parametric study are presented to clearly define the key design variables and their associated 
influence. In general the important parameters were identified as misfit, mean radius, thickness, material property 
and/or load gradation, and speed; all of which must be simultaneously optimized to achieve the “best” and most 
reliable design. Also, the important issue of defining proper performance/merit indices (based on the specific stored 
energy), in the presence of multiaxiality and material anisotropy is addressed. These merit indices are then utilized 
to discuss the difference between flywheels made from PMC and TMC materials with either an annular or solid 
geometry. 

Finally two major aspects of failure analysis, that is the static and cyclic limit (burst) speeds are addressed. In 
the case of static limit loads, upper, lower, and out-of-plane bounds for disks with constant thickness are presented 
for both the case of internal pressure loading (as one would see in a hydroburst test) and pure rotation (as in the case 
of a free spinning disk). The results (interaction diagrams) are displayed graphically in designer friendly format. For 
the case of fatigue, a representative fatigue/life master curve is illustrated in which the normalized limit speed versus 
number of applied cycles is given for a cladded TMC disk application. 


SYMBOLS 


Material Parameters 


h 

e t 

K 


Or 


P 

v, 


a ,a 

U C U T 


longitudinal Young’s modulus 
transverse Young’s modulus 
transverse yield stress under shear loading 
longitudinal thermal expansion 
transverse thermal expansion 


ratio of anisotropy, i.e., I ^ 


longitudinal Poisson’s ratio, i.e., (e r /S q, load is in 0 direction) 


ratio of longitudinal to transverse ultimate stresses, i.e.. 


\® u t J 


material density 

ultimate tensile strength (UTS) in the longitudinal and transverse directions respectively. 


g y l longitudinal yield stress 
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Stresses and Strains 


Dy second order direction tensor 

d i unit vector denoting local fiber direction 

F ( ) limit function (static fracture, endurance limit and normalization surfaces) 

Sy deviatoric stress tensor 

Gy Cauchy stress tensor 

G r Zj. radial stress and strain components respectively 
G 0 , 80 tangential stress and strain components respectively 

Subscripts and Superscripts 


( ) L X ) T longitudinal and transverse properties (static fracture, endurance limit and normalization surfaces) 
respectively 

(•) denotes differential with respect to time 


Geometry 


a 

b 

h 

R m 

R i 

r 

u 

5 

X 


inner radius 
outer radius 
height 
mean radius 

radial location of the ith interface 
radius 

radial deflection 

misfit between concentric disks 

width of a given disk 


Forces 


P {n internal applied pressure 

P QU[ external applied pressure 

p Generalized Pressure Loading: 

(a) When internal pressure loading applied; equal to internal pressure P m 

(b) When rotational loading is applied; equal to the equivalent centripetal pressure load, i.e., 

P= pco 2 ( b 2 - a 2 )/2 

T temperature at the inner diameter, r = a 

T h temperature at the outer diameter, r = b 

AT radial temperature distribution 

(0 speed of rotation 


1.0 INTRODUCTION 
1.1 General 

Over the years, there has been an increasing demand for effective energy storing systems (ref. 1). To meet such 
needs, many alternative systems have been proposed, as listed in table I, with composite flywheel energy storage 
systems rising to the top of the list of candidates, because of there high power and energy densities with no fall-off 
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in capacity under repeated charge-discharge cycles. In recent years, a host (e.g., refs. 1 to 3) of different flywheel 
system designs (composed of various fiber-reinforced material systems, i.e., metallic composites, polymeric com- 
posites, etc.) have been identified with the goal of improving the performance of these systems. These flywheel sys- 
tems can be categorized into principally three design topologies: (1) preloaded, (2) growth matching, and (3) mass 
loaded, with each approach having its own set of advantages and disadvantages. 

Energy storage in all these flywheel designs is basically sustained by a spinning mass, called a rotor. This 
important “rotor” component has therefore been the subject of research investigations over many years, generating 
extensive literature on its many analysis and design related aspects. At different stages of developments, a number of 
comprehensive state-of-art reviews (refs. 1 and 2) have been completed. For instance, representatives of these in the 
area of material-system selections for flywheel devices are given in Genta (ref. 1) and DeTeresa (ref. 4) whereas 
optimization for integrated system design are given by Bitterly (ref. 3). 

In the following section, a brief summary is given for the state-of-art that is presently available on the analysis- 
related aspects of the rotor-type components (also disks, multiple-disks, etc.). However, we emphasize at the outset 
that the present review is not intended to be exhaustive; rather we are primarily concerned with works and studies on 
the important factors that should be considered in any valid assessment of the performance (both at working/service 
and at limiting operational conditions) of flywheel systems. To this end, we will group the pertinent list of refer- 
ences in tables II to IV according to three major considerations; i.e., (1) factors related to additional loading condi- 
tions beyond those of rotation (e.g., internal pressure, temperature, and interference fits, etc.), (2) factors related to 
material behavior (e.g., isotropic versus anisotropic, elastic versus inelastic, etc.) and (3) factors related to geometry 
(e.g., uniform versus variable thickness, solid versus annular). This background then puts into perspective the subse- 
quent parametric studies and related conclusions in the later part of the paper. 


1.2 Background Literature 

Starting as early as 1906 the rotating isotropic disk has been studied by Grubler (ref. 6) followed by Donatch 
(ref. 7) in 1912. One of the first “modem” dissertations in analyzing the flywheel rotor was the seminal work done 
by Stodola (ref. 8), whose first translation to English was made in 1917. In 1947 Manson (ref. 9) published a finite 
difference method to calculate elliptic stresses in gas-turbine disks that could account for the point-to-point varia- 
tions of disk thickness, temperature, and material properties. Millenson and Manson (ref. 10) extended the method to 
include plastic flow and creep in 1948. Both methods have been widely used and extended by industry (ref. 11). 
More recently, Genta (ref. 11) modified Manson’ s method to include orthotropic materials. The 1960s and 1970s 
lead to numerous other serious efforts in analyzing the rotor, and introducing different designs for the flywheel, with 
the onset of composite material development giving added impetus. A detailed review of the rotating disk problem 
up through the late 1960’s is given by Seireg (ref. 12), whereas, Habercom (refs. 13 to 15) provides a collection of 
abstracts and reference of government-sponsored flywheel related research (late 1960's until the late 1970’s). The 
past 40 years showed a tremendous amount of work done concerning the flywheel, e.g., see Genta’s (ref. 11) exten- 
sive review, including a historical perspective, on research work conducted up to the early 1980’s related specifi- 
cally to flywheel systems. 

For convenience in presentation, we have grouped in tables II and III, respectively, the different representative 
works with regard to the two major material behavior classifications treated; i.e., elastic-isotropic and elastic aniso- 
tropic. These classifications in fact have had a long-standing history in the available literature. The more complex 
nonlinear counterpart (i.e., plasticity, creep, fracture, etc.) is of more recent origin, and consequently comparatively 
smaller in size, and primarily focused on isotropic material behavior. A sampling of work on these latter nonlinear 
areas is given in table IV. 

With reference to tables II and III we note that most of the listed investigations have focused on the disk of uni- 
form thickness under rotation. Alternatively, some investigators such as Kaftanoglu (ref. 5), Abir (ref. 41), Hunting- 
ton (ref. 70), Nimmer (ref. 90), Daniel (ref. 94) and Portnov (ref. Ill), considered the effect of misfit in addition to 
disk rotation as they studied the multirim disk case. Others like Leopold (ref. 18), Potemkina (ref. 31), Prasek 
(ref. 39), Shanbhag [46], Chakrabarti [61], Gurushankar [66], Dick [74], Baer [82], Misra [100], Ferrero [102] and 
Portnov (ref. 108), considered the temperature effect as well as the disk rotation case. 

Considering the list in table IV, we refer to Shevchenko (ref. 1 19), Weber (ref. 120), Pisarenko (ref. 121), Reid 
(ref. 122), Gururaja (ref. 127) and Gueven (ref. 136), for a representative nonlinear elastoplastic analysis of rotating 
disks, and to Lenard (ref. 147), Gupta (ref. 150), and Mukhopadhyaya (ref. 152), for a number of time-dependent 
(viscoelastic) solutions. 
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1.3 Objectives and Outline 


In summary, tables II to IV with their broad classification of previous work, convey the complexity of the prob- 
lem, in that a very large number of conditions and factors must be considered. In particular, each of the listed inves- 
tigations have necessarily focused on some specified conditions, particular material behavior and related 
conclusions. An urgent need therefore exists for a compilation, in as concise a form as possible, of the extensive 
results, major trends observed and the many conclusions reached in these previous studies. To the authors knowl- 
edge a study of this type is presently lacking; in which the problem in its “totality’' has been considered. As a first 
step toward this end, an annotative review study is conducted herein, with the major restriction being that of disks 
with constant thickness. In particular, emphasis is on combined loading (rotation, pressure, misfit, and temperature) 
and the effect of material anisotropy (with its great impact on stress distributions and any ensuing singularities in 
stress, strength and failure modes, etc.). To facilitate the presentation of the well-established elastic anisotropic 
solution procedure (single and multiple disks) specific reference to the numerous previous works cited in tables II 
and III will not be individually mentioned. However, we emphasis at the outset that all the conclusions made herein 
are in agreement with these previous works, except of course the notation presently employed. In addition to this 
review, our second primary objective is to present a number of new analytical and numerical results with regard to 
limiting conditions of burst (pressure as well as rotation), fatigue, etc. These will then enable us to reach a number of 
important conclusions in parallel under both service conditions (stress analysis), as well as at the overload conditions 
for “lifing” studies of the problem. 

In an outline form, the remainder of the paper is described as follows. In section II, we give a detailed deriva- 
tion of the analytical model used in the elastic stress analysis for single/multiple, annular/solid, anisotropic/isotropic 
disk systems, subjected to both pressure surface tractions, body forces (in the form of temperature-changes and rota- 
tion fields) and interfacial misfits. Section III contains the results of extensive parametric studies and important 
observations obtained therefrom. In Section IV, we address the important issue of defining proper performance/ 
merit indices (based on the specific stored energy), in the presence of anisotropy. Section V deals with two major 
aspects of failure analysis, i.e., static and cyclic, providing the corresponding limit (burst) load and representative 
fatigue/life limit speed curve, respectively. 


2.0 ELASTIC STRESS ANALYSIS VIA ANALYTICAL MODEL 

2.1 Single Disk Solution 

2.1.1 Annular Disk . — Here, as classically done the problem of a quasi-static rotating disk subjected to both ther- 
mal and mechanical boundary conditions is summarized by stating the governing equations required to define a 
well-posed boundary value problem. We begin with the equation of equilibrium: 


da r | ~ c e 

dr r 


+ pco 2 r = 0 


( 1 ) 


or 


<Jg =^-(ra r ) + pco 2 r 2 (2) 

dr 

in which the centrifugal “internal force” (pco 2 r) is included as a body force term. Next the pertinent compatibility 
equation, 


£r 



is obtained from the well known strain — displacement conditions 


( 3 ) 
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and 


du 

dr 


( 4 ) 


u 

£d=~ 

r 


(5) 


so that only the specific constitutive equations to be utilized, need be defined. Given that the problem is an 
axisymmetric, transversely isotropic disk (with the fiber directed along the circumferential direction of the disk 
fig. 1), the plane stress constitutive relations are as follows: 


c r 

£,■ = — 


v L°e 


+ cCf AT 


( 6 ) 


s e = - Vl<j + a L AT 


(7) 


where AT is the assumed linear temperature gradient across the disk and is defined as 

T b -T n 


AT = T„-T n +- 


b-a 


— ( r-a ) 


( 8 ) 


where a and b define the inner and outer radius of the disk, T Q and T b are the imposed temperature at the inner and 
outer radius respectively, and T 0 is the reference temperature of the disk. 

Substituting the constitutive relations equations (6) and (7) into the compatibility expression, equation (3), and 
utilizing equilibrium, equation (2); the following governing equation for the radial stress component is obtained: 


where 


2 dr Gy , ^ dc r 


dr 


+ 3r — + (1 — — )<7 — Q(r) 
dr p 


(9) 


P = ^r- and £?('*) = T~~ 
E l b-a 


r ^L^a ~T b )(2a L -a T ) 

+a^E L {T b - T 0 ){ a i -a r ) + r“pco“(3 + v L )| 

-b^E L (T a -T Q )(a L -a r ) + /* 2 p(o 2 (3 + v L )} 


(10) 


Implied in equation (9) is the assumption that the material parameters (i.e., E L , E T , v L , v r , and G L , a L and Oj) are 
not strongly temperature dependent and therefore can be taken to be independent of radial location. Inclusion of this 
temperature dependence would only contribute higher order effects in the solution and distract attention from the 
primary emphasis of examining the effects of rotation. Also, the only practical values for p are those less than or 
equal to one, as the fiber stiffness is always greater than the matrix (thereby resulting in E L > E T ). 

The general solution of this Euler-Cauchy governing equation is: 


c J / .=C 1 ;* mi +C 2 /^ 2 +G 1 (r)-e 2 (r) 


(ID 


where 
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m x 

Ql(r) = f r~ (m ' +l> Q(r)dr 

m\ - m 2 J 

m 2 r 

Q 2 (r) = f r~ (mi+l) Q(r)dr (12) 

mi - m 2 J 

mi = -1 + X, m 2 = -1 - X and Z = 



and C { and C 2 are constants to be determined from the applied mechanical boundary conditions. 

The general solution for the tangential stress is obtained by merely substituting equation (11) back into (2): 

o e =Ci(m { +l)/* mi +C 2 (m 2 +l)r m2 +(m { + l)Q { (r) - (m 2 + l)Q 2 (r) + pco 2 r 2 (13) 

Whereupon, the radial deflection is derived by substituting equations (11) and (13) into (7) and the resulting expres- 
sion into equation (5): 


u r 



VL-°r 

E L 


+ oc L 



T b -Tg 

b-a 


C r-a ) 


) 


(14) 


To obtain the constants C x and C 2 appearing in equations (11) and (13), we impose the following boundary 
conditions along the inner and outer radius of the disk, i.e., 


«y r )r=a = -Pin 
(&r \=b ~ ~Pout 

where P in and P QUt are the internal and external pressures, respectively. Solving the resulting 2x2 matrix equation 
results in the following expressions: 


Q = 



1+ 

[Pin +C?2( a )-C?l( a ))- fc 





(15) 


and 


p.'iP 
a ' ' b ’ K 


C,=- 


-ab W (P in + Q 2 (a) - Qt (a)) + bJ P (P out + Q 2 ( b ) - Q, (b)) 


7 II 


a » 


(16) 


where we may separate the thermal and rotational terms such that 

Ql(r) = Qf(r) + Qi R (r) and Q 2 (r) = Q% (r) + Q? (r) (17) 

with the generalized thermal forces being defined as, 
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( 18 ) 


Cf(D— 1 


2 b 


F tn n \ a ( T b ~ T o'> + b(.-T a + T 0 ) rE L (T a -T b )( 2 a L -a T ) 
t L (a L -a T ) p + j prr 

- ,+ ^H M r b) 


Ql(r) = 


1 




^L^ a L ~ a T ) 


a(T b -T 0 ) + b(-T a +T 0 ) rE L (T a -T b )(2cc L -a T ) 


-1-J^ |(«-^) 


• + 


~ 2 ~-lj K^-b) 


and those due to rotation as. 


fiftr). ^ 2(3+V A gf(r)= ^ 2(3 + V ^> 


2 |I 

Ip 


3 + p 


2 


3 P 


(19) 


2.1.2 Solid Disk . — For the special case of a solid disk, the radial stress must remain finite at the center of the 
disk, therefore given either an isotropic or anisotropic material, it can be immediately seen that C, in equation (11) 
must be equal to zero (thus both displacement and stress fields at the center will be finite), such that, 

O r =C x +Q { {r)-Q 2 {r) (20) 


for the isotropic case, and 


<J r =C l r m ' +Q l (r)-Q 2 (r) (21) 

for the anisotropic case. Now assuming that an applied external pressure P QU[ exist, the constant C l is determined to 
be: 


Q ~ P out + Ql W “ Q\ W 


( 22 ) 


for the isotropic case, or 


Ci=^[P out +Q2W-Qm (23) 

for the anisotropic case. Where in the definitions for Q\, Q r and Q\, Q 2 given above in equations (18) and (19), a 
(the inner radius) is taken to be zero. Similarly, the tangential stress is 

Gq = C\ + Q\ (*0 + Qi (r) + pco 2 r 2 (24) 


for the isotropic case, and 


o e = (w, + 1 )C 1 /- m ' + (m, + 1)2, (r) - (w 2 + l)g 2 (r) + pco 2 r 2 


(25) 
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for the anisotropic case. Note, for the isotropic case, Q x and Q 2 , can be greatly simplified by taking (3=1, with the 
resulting expressions being given in equation (29), see section 2.4.1. 


2.2 Multiple Number of Concentric Disks 

Extension of the formulation to “/?” number of annular concentric disks, with possible misfit (see fig. 2), is sim- 
ply accomplished through application of appropriate boundary conditions at the interface of each disks, so as to 
ensure continuity of radial stresses and displacements. Continuity of radial stress demands the following: 

cr, .(/• = /}) = />, 

a l] (r = r 2 ) = a, i (r = r 2 )= P 2 

a r 2 ( r = r 3 ) = °V 3 (r = r 2 ) = P 2 

. “ (2o) 

0=r„ ) = <7 rn (r = r „) = P n 
<J r„( r = r n+0= P n+l 

where P j and are the applied internal and external pressures, respectively; whereas the remaining continuity 
conditions in (26) will then lead to the determination of the remaining interfacial pressures P 2 , P y .. P n (which are 
redundant unknowns, for the common case of imposed (known) interference fit). Similarly, kinematic constraints are 
introduce at the interfaces corresponding to the imposed misfit (5) at each interface. These are: 

u r 2 ( r = r 2> ' u r, ( r = r 2> = 5 1 

u r 3 < r = r 3> - u r 2 ( r = r 3> = S 2 

u r 4 (r = r 4 ) - u rj (r = r 4 ) = <5 3 (27) 

u r„ (r = r n > - u r n .| ( r = r n> = ^n-l 

Where 5j, 8 0 , ...8^ are the imposed misfits between disks that will in turn generate interface pressures. Combining 
equations (26) and (27) with equations (11) and (14) specifies a system of 2 n equations with 2 n number of unknown 
constants; i.e.: 


[*]{<:}={*} 


(28) 


where 

{ C } is the vector of unknown constants; {C} = jc* , c\ , c\ , cf , • • • , C" , C£ J 

[AT] is the geometry and material matrix 

[R] is the vector of applied forces and/or misfit parameters 


2.3 Verification of the Solution 

As indicated in the introduction, the analysis of isotropic and anisotropic disks subjected to rotation, internal 
and external pressure, temperature variation and misfit are not new (refs. 1 to 4). Consequently, the present unified 
solution can be, and has been, validated through comparison with a variety of other solutions, see table V for a brief 
summary. Clearly all aspects of the solution have been verified with at least one other independent numerical or 
analytical solution. In all cases the present solution agrees identically with all other analytical solutions and for all 
practical purposes identically with all FEA results as well. Figures 3 to 5 illustrate three of the many comparisons 
made between the present formulation results and other previous analytical results. In figure 3 the thermal and inter- 
nal/extemal pressure analysis is compared against a previous analytical solution derived by Arnold (ref. 175), with 
the agreement being excellent. In figure 4 the analysis of an anisotropic disk subjected to pure rotation is compared 
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with Lekhnitskii’s solution (ref. 58), again the agreement is excellent. Another example demonstrates the ability to 
include misfit within the formulation, see figure 5. Here comparison with Ugural and Fenster’s (ref. 176) solution of 
two isotropic disks subjected to a misfit of 0.1 mm is made. The calculated contact pressure is, 12.305 MPa, while 
the tangential stresses in the outer and inner cylinder are: 

G 0 (/* = 200 mm) = 56.055 MPa 
G 0 (r = 250 mm) = 43.750 MPa 


and 


G 0 (r = 150 mm) = -56.250 MPa 
G 0 (r = 200 mm) = -43.945 MPa 

respectively. These values agree exactly with those calculated by Ugural and Fenster (ref. 176) 


2.4 Observations Regarding The Analytical Solution 


Before conducting a parametric study to identify the important parameters in the design of flywheel systems, let 
us make a number of observations regarding the current analytical approach for conducting stress analysis by exam- 
ining a number of special cases. 

2.4,1 Material Isotropy . — For the special case of material isotropy (i.e., (3=1), the character of the governing 
equation (see eq. (9)) will change as the zeroth order term will disappear and consequently, the particular solution, 
i.e., equation (12), will be appropriately modified. The resulting expressions for Q { and Q 2 coming from 
equation (12) and including both thermal and rotational loading are: 


Ql(r) = r 

/ 

Q 2 (r) = r 


E(T a -T b )a 
2 (a-b) 


E(T a -T b )g 
~6(a - b) 


rpco~( 3 + v) 
-4 


rpco 2 ( 3 + v) > 
-B , 


(29) 


2.4.1 Thermal Isotropy . — Here let us consider the case when only thermal loading is applied to the disk. Under 
this condition it is apparent from equations (10) to (19) that, in general, the stress distribution (tangential and radial) 
is dependent upon geometry, material properties (i.e., strength of anisotropy), mismatch in thermal expansion coeffi- 
cients, and the magnitude of the thermal load itself. In the case of thermal isotropy (i.e., when the tangential (a L ) 
and radial (ct^) thermal coefficients are equal) the stress distribution is primarily dependent only on (3 (the ratio of 
moduli) and E L , as the constants C { and C 2 become a function of Q 2 - Q v that is: 


r n rE L (T a -T b )cc 

V 2 ~vi ~ r~ 

(a-b)( 4--) 


(30) 


Consequently, under a uniform temperature distribution no stresses are developed as one might expect. Also in gen- 
eral, from equations (10) and (14), it is apparent that the longitudinal thermal coefficient is the dominant parameter 
with the longitudinal modulus dictating the magnitude of thermal stresses being developed. 

2.4.2 Numerical Singularities . — It is well know that nonphysical numerical singularities can exist within an 
analytical solution, and these vary depending upon the approach taken. In the above derivation the final governing 
equation to be solved is stress based, therefore we will examine the resulting radial stress under various special con- 
ditions so as to identify all potential singularities. As the solution is limited to elastic behavior the concept of super- 
position applies; consequently individual-loading conditions can be examined separately. This will greatly simplify 
the current task. The following four possible loading cases (rotation, pressure, thermal and interfacial misfit) were 
taken under consideration. 
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Rotation Only : Given a single, annular, disk subject to only rotation, the radial stress equation (1 1) becomes: 



Clearly, a singularity would exist in the solution whenever p approaches 1/9 (i.e., whenever E L = 9 E T ). This type of 
singularity has been found before in other solutions, such as Lekhnitskii (ref. 58) and more recently for the case of 
solid disks, as discussed in reference 177, for values of p > 1. Once again confirming the accuracy of the present 
formulation. 

Pressure Only : Solving for the single, annular, disk case subjected to both inner and outer pressure {P- and 
P ouV respectively), the following equation for the radial stress can be derived: 
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(32) 


For any finite disk, a can not equal b in equation (32), and consequently there is no possible singularity for this type 
of loading. 

Thermal Loading Only : Again considering a single, annular disk, subjected to a pure thermal loading condition, 
the derived expression for the radial stress becomes somewhat lengthy; consequently, for convenience, only the 
denominator (which would be the source of any numerical problem) of the troublesome term will be shown here. 
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Clearly the solution has two singularities, namely whenever p = 1/4 (i.e., E L = 4 E T ) or p = 1 (i.e., E L = E T , which is 
the isotropic case). 

Misfit only : The simplest problem containing a misfit is that of two disks with a misfit, 8, imposed between the 
inner and outer disks, see figure 6, where a and b are the inner and outer radii, respectively and c designates the 
radial location of the misfit. Solving for the radial stress in both the inner and outer disks, respectively, we get the 
following expressions: 
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Inner disk 
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Outer disk 



( 34 ) 


(35) 


It is obvious from examining equations (34) and (35) that for the case of misfit only, no singularities exist over the 
domain. 

Given the results of the above four loading cases it is apparent that a pattern has emerged, thus prompting a 
generalization regarding the potential for the existence of numerical singularities. One can expect singularities to be 
present whenever body forces are imposed upon a given solution : however, when the excitation is arising due to the 
application of imposed boundary conditions (i.e., applied tractions or displacements) no singularities should be 
expected to occur. 


3.0 INFLUENCE OF KEY DESIGN PARAMETERS 

Having developed and verified a general analytical solution for multiple, anisotropic, concentric disks, with 
misfit, subjected to rotation, pressure and thermal loading, we are now in a position to identify and examine the 
influence of important design parameters such as: mean radius, thickness, the amount of misfit, rotational speed, and 
material gradation. Two classes of anisotropic (composite) materials are primarily employed throughout this study. 
The first class, polymeric matrix composite (PMCs) systems, is represented by a 60 vol % fraction graphite/epoxy 
system. The second class, titanium matrix composite (TMCs) systems, is represented by a 35 vol % fraction SiC/ 

Ti- 15-3. Subsequent to examining a single disk composed of a uniform volume fraction of material, examination of 
material gradation and misfit effects will be undertaken by considering three concentric disks with varying fiber 
volume fractions. 

Flywheels composed of PMC material systems are currently being investigated under the NASA funded Rotor- 
Safe-Life Program (RSL); a collaborative govemment/industry program designed to establish a certification proce- 
dure for composite flywheels. Consequently, this material class will dominate the parametric study. The associated 
material properties for the two classes of composite systems (and their individual fiber and matrix constituents) used 
in this study are listed in tables VI and VII, where the composite properties were obtained using the micromechanics 
analysis code based on the generalized method of cell, MAC/GMC (ref. 178). All properties given are at room tem- 
perature and assumed to be independent of temperature. The flywheels to be considered in the RSL program are 
expected to have a mean radius of approximately 4 in. and a thickness of approximately 3 in.; therefore these geo- 
metric properties will constitute our baseline case. 
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3.1 Single Disk Case 


Considering the case of a single annular disk, subjected to pure rotation, the influence of the three key material 
and geometric parameters (i.e., the degree of material anisotropy ((3), mean radius (R m ) and thickness (p/.A,)) on the 
radial and tangential stress distribution will be investigated parametrically. Analytically, these parameters can be 
clearly seen by normalizing both the radial (see eq. (31)) and tangential stress (see eq. (2)) with respect to the den- 
sity, rotational speed, and outer radius, b , that is: 





(36) 



V / 

where the inner and outer radii are defined as 


b = R m +\ 


and 


a 
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2 R 

~ JV /77 

1 + — 
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so that g(p, R m , A,, v L ) and/(p, R m , X , v L ) will be denoted as the radial and tangential distribution factors, respec- 
tively, for an anisotropic annular disk. 

Figures 7 and 8 depict these factors as a function of normalized radial location, r/R m , for two limiting geometric 
representations, i.e., a thick disk (X/R m ~ 2.0) and a thin disk (X/R m « 0.125), given a slightly anisotropic (p = 0.7) 
TMC, and strongly anisotropic (P = 0.05) PMC material description, respectively. Further, in table VIII the maxi- 
mum values of the modified normalized radial and tangential stress are given for a range of thickness and mean 
radii. From these results one can draw a number of conclusions. 
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The first being that both the radial and tangential stresses increase as R m and X are increased (see table VIII). 
The second is that the degree of material anisotropy (p) impacts the magnitude of increase in the radial stress com- 
ponent more so than the tangential stress component, and this influence increases as the thickness is increased; 
whereas, it diminishes as the mean radius is increased. In fact, once the disk is sufficiently “thin,” i.e., X/R m < 0.125, 
the degree of material anisotropy becomes irrelevant. Thirdly, the location of the maximum stress is also greatly 
impacted by the degree of anisotropy. For example in figure 7, we see that for the strongly anisotropic PMC material 
system the maximum tangential stress location varies greatly depending upon mean radius; whereas for the slightly 
anisotropic TMC material system, see figure 8, the maximum tangential stress is always at the inner radius. Simi- 
larly, the location of the maximum radial stress shifts from approximately 0.75 R m to R m as the mean radius 
increases for a TMC (see fig. 8). Whereas for a PMC (see fig. 7) material the maximum location moves from 
1.25 R m , for small mean radii, towards the center (i.e., R m ) as the mean radius, R m , is increased. These observations 
were substantiated by the work of Gabrys and Bakis (ref. 179) wherein they designed and tested a flywheel that 
employed a much softer urethane matrix, which precluded the development of a large radial stress. This then 
allowed the use of thick flywheels, wherein the maximum hoop stress was shifted so that it occur near the OD; 
thereby theoretically leading to a less catastrophic failure mode and fail-safe flywheel design. Fourthly, irrespective 
of the degree of material anisotropy the most efficient use of material is made when one considers the case of rela- 
tively thin disks rather than thick ones, see figures 7 and 8. This will be discussed further in the subsequent material 
and configuration selection section. Finally, it is obvious from equations (36) and (37) that the density (p) of the 
material impacts the actual magnitude of the stress field (both radial and tangential components) linearly, while the 
rotational speed (co) impacts both components quadratically. 

The radial and tangential stress distributions as a function of normalized radial location for pure thermal loading 
and internal pressure loading are shown in figures 9 and 10; assuming the baseline material parameters correspond- 
ing to the 60 vol % fraction PMC system in table VI for a relatively thick disk (R m = 4 in. and A,= 3 in.). From 
figure 9, we see that temperature (whether uniform or gradient) induces compressive hoop stress in the outer portion 
of disk and tensile hoop stress in the inner portion. However, in the presence of a linearly varying temperature field, 
the slope (positive or negative) of the gradient of temperature will shift the maximum location of the tensile radial 
stress from one side of the disk to the other, as well as change the curvature of the hoop stress distribution and the 
magnitude of the maximum tensile hoop stress. 

In figure 10 the influence of applying pressure along the inner surface of the disk is illustrated. This type of 
loading distribution would be consistent with the presence of a solid hub (shaft), albeit in practice it is clear that the 
magnitude of this internal pressure would be a function of the density, size and speed of rotation of this solid hub. 
Consequently, in figure 10 both the radial and tangential stress profiles are normalized with respect to the applied 
internal pressure. The basic influence of internal pressure is to induce a compressive radial stress state and tensile 
hoop stress throughout the annular disk, with a maximum in both stress components occurring at the inner radius. 
Further, discussions regarding the influence of intemal/extemal pressure and temperature profiles for single cylindri- 
cal or disk configurations can be found in reference 175. 

Now the superposition of these three classes of loads (thermal, internal pressure, and rotation) is illustrated in 
figures 1 1 and 12; where both radial and hoop distributions are shown for a relatively thick disk (^ = 3, R m = 4) 
rotating at 15 000 rpm (fig. 11) and 60 000 rpm (fig. 12) subjected to both uniform and gradient temperature profiles 
and an internal pressure of 1 ksi. Note how for the lower speed case (fig. 11), the overall conclusions for the com- 
bined loading situation are similar to those of the pure thermal cases; with the exception of tensile hoop stress being 
present even in the outer region of the disk for both the uniform and inner gradient case. Tensile hoop stresses are 
present in these two cases because the tensile hoop stress generated, due to rotation, are sufficiently high to over- 
come the thermally induced compressive hoop stresses in the outer region. In addition we see immediately that at the 
inner radius the radial stress becomes compressive due to the applied internal pressure, while it is zero at the outer 
surface, as it must be. Figure 12 clearly shows that as the rotational speed is increased to 60 000 rpm the stresses due 
to rotation dominate those induced by either the internally applied pressure or the imposed thermal profiles. Note, 
however, that the hoop stress distribution becomes almost uniform across the disk, thus making more efficient use of 
the material. 
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3.2 Three Concentric Disks 


In the RSL program the flywheels of interest are being manufactured with a preload design in mind, that is the 
flywheel is composed of a number of concentric rings press-fit together with a specified misfit between each ring. 
Therefore, the resulting laminated flywheel has a compressive radial stress at the interface of each concentric ring. 
Consequently, the desire here is to investigate the influence of a key design parameter (i.e., misfit, 5 between disks) 
on the corresponding radial and tangential stress distribution within the laminated system. This parameter is in addi- 
tion to the two previous geometric parameters: R m - mean radius and X - ring thickness in the radial direction. When 
examining the influence of 8, tracking of the radial stress distribution is important, since a positive value at a given 
interface will constitute (by definition) disk separation. Similarly, if the radial stress exceeds some critical value, 
delamination within a given disk may also take place. Alternatively, the tangential (or hoop) stress is significant in 
that it will dictate the burst failure of a given ring (disk) within the system, which may or may not (depending upon 
the location) lead to an immediate catastrophic failure of the entire system. The limit speed of a given, 
circumferentially reinforced, disk is dictated by the circumferential (or hoop) stress, thus excessive hoop stress will 
result in fiber breakage and a subsequent total loss in load carry ability. Consequently, this failure criterion is 
defined as the burst criterion. 

3.2.1 Uniform Material Properties . — To simplify the study, only three (each with the 60 vol % fraction PMC 
system given in table VI) concentric disks will be investigated with a baseline geometry similar to the previous 
single disk case, that is, an R m = 4 and a total thickness (i.e., X^ + X~> + X 3 ) equal to 3 in. A three disk minimum is 
required so that the influence of two misfits can be examined, that is the misfit between disk 1 and disk 2, 8j, and 
that between disk 2 and disk 3, 8 2 , see figure 13. When considering applying a misfit at an interface one could 
specify a similar misfit directly for each interface or alternatively specify a given mismatch in hoop strain for all 
interfaces. The first option would effectively result in applying different amounts of hoop strain at each interface, 
whereas, the second option would result in different misfits being specified for each interface. In an attempt to iso- 
late the interaction of mean radius and misfit, the constant applied hoop strain approach to selecting the appropriate 
misfit at each interface is taken. Consequently, each displacement misfit applied is proportional to the radial location 
of the corresponding interface (i.e., 5, = erfl.) such that a constant hoop strain (£ /? ) at each interface is imposed. In 
this way the impact of the misfit will be consistently felt even when the mean radius is increased. 

Assuming only rotational loading with a speed (CD) of 30 000 rpm, analysis results for the disk system are shown 
in figures 14 and 15, when the overall mean radius and misfit hoop strain are varied as follows: R m = 2, 4, and 6 in. 
and z h = 0., 0.1 and 0.2 percent. In figures 14 and 15, the radial stress versus radial location and tangential stress 
versus radial location are shown respectively, for each of the above cases. Examining figure 14 it is apparent that in 
the case of no misfit (e /7 = 0, see fig. 14(a)) all radial stresses are tensile. Recalling from figures 1 1 and 12, it is clear 
that even if internal pressure was applied to offset these rotationally induced stresses only the inner most surface 
would experience compressive radial stress, while the remainder of the disk would remain in a state of tension. The 
extent of the compressive region is dependent upon (as shown in fig. 10) the magnitude of the imposed inner pres- 
sure. However, increasing this inner pressure (e.g., in the case of a mass loaded design) is of only limited value as 
this compressive region rapidly becomes tensile as one proceeds toward the outer diameter of the disk system. 

Alternatively, if one imposes a prestress through application of a misfit at the interfaces between disks, a com- 
pressive radial stress state can be induced throughout. These interfacial compressive radial stresses (denoted in fig- 
ures 14(b) and (c) by an /j and / 2 ) then increase as the misfit hoop strain is increased for a given mean radius. Note 
that even though the imposed hoop strain mismatch is the same at all interfaces, within the system, the first inter- 
face, 7p is more compressive as it benefits from the compressive state at / 2 as well. Furthermore, figure 14 clearly 
demonstrates that as the overall mean radius of the system is increased, the effectiveness of the imposed hoop strain 
is reduced; even to the extent of not inducing any compressive stresses. 

Figure 15 illustrates the tangential stress distribution for various mean radii and magnitudes of hoop strain mis- 
match (misfit) between the disks. Examining figure 15, it is apparent that, (1) the location of the maximum tensile 
hoop stress (and thus the initiation of fiber breakage and failure) changes and (2) the magnitude of the hoop stress 
increases as the mean radius is increased, given a specified amount of hoop strain mismatch. In addition, in the pres- 
ence of misfit, the tangential (hoop) stress profile is discontinuous and forms a stair step pattern, with each subse- 
quent disk within a given system, typically having a higher stress state than the previous. Consequently, given an 
interfacial misfit and the given material anisotropy, one can be assured that failure would initiate at the inner radius 
of the outer most disk, provided the overall system is not relatively thin, i.e., XJR m > 0.3. Nevertheless, it is apparent 
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that a complex relationship between the various geometric parameters, 8, X, and R m exists (as suggested by eqs. (34) 
and (35) for the case of two disks). 

Selecting a proper disk thickness configuration is another important aspect that was investigated. Figure 16 
shows three arrangements with different individual disk thickness but all adding up to the same overall total thick- 
ness. Examining the figure it can be concluded that the most beneficial arrangement is that of two thinner inner rings 
surrounded by a thicker outer ring. As this arrangement provides a maximum compressive radial stress state while 
giving a lower maximum hoop stress within the outer ring. Clearly, the thicker outer ring could fail (as could all 
other configurations) by delamination due to excessive tensile radial stress. However, for this particular case delami- 
nation and burst would be confined to the outer ring, which should be a more fail-safe design than either of the other 
two configurations examined. 

3.2.2 Nonuniform Material Properties and/or Loads . — The goal of any design is to maximize the allowable rota- 
tional speed before either a disk separation (in a preload design) or burst criterion is exceeded. Although, in equa- 
tions (34) and (35), both disks were assumed to be comprised of the same transversely isotropic material, it is 
apparent that not only is the magnitude of misfit important in the determination of the radial stress at the interface 
but so too are the material properties. Consequently, this led us to investigate the influence of gradation of the mate- 
rial properties within the laminated system. As a result, the six cases identified in table IX were examined. The first 
three cases being PMC (strongly anisotropic) laminated systems and the last three TMC (slightly anisotropic) sys- 
tems. An additional motivation (besides the degree of anisotropy) behind investigating an idealized TMC system as 
well as a PMC system, is the ability to examine the trade-offs between density and stiffness/strength. Since in a 
PMC system, as the volume fraction of reinforcement is increased the density (weight) of the material as well as 
stiffness and strength also increase. Alternatively in a TMC as the fiber volume fraction is increased the density of 
the material is reduced, while the stiffness and strength are increased. Albeit, this later density consideration will be 
small and in practicality higher order when compared to the variation in stiffness and strength. Even so, the induced 
body forces due to rotation, and thus influence of gradation effects is different between the two classes of materials. 

Case one is identical to the baseline material examined thus far; consequently, case two through six will be sub- 
jected to the same conditions. These are a rotational speed (co) of 30 000 rpm, variation of mean radius (R m = 2, 4, 
and 6 in.) and imposition of constant hoop strain (e h ) at each interface of 0, 0.1 and 0.2 percent (i.e., variable misfit, 
8, per interface such that 8, = E h Rj). Again the representative model TMC material properties are given in table VII. 
The results of these analyses are presented in figures 17 to 20, where the radial stress versus radial location and tan- 
gential stress versus radial location are shown, respectively, for all three cases. 

Examining these figures it is clear that material and density gradation play a major role in the resulting radial 
and tangential stress distributions; the geometric and strength of anisotropy effects being similar in nature to those 
discussed at length previously. In particular, note how the radial compressive stress at the disk interfaces can be 
increased or completely eliminated depending upon the imposed material gradient. Similarly, the magnitude of the 
stair steps in the tangential stress as well as the maximum value and location of this stress are also greatly impacted 
by gradation. For example, note the difference between case two and three (compare, figs. 18(b) and (c)) in that for 
case two burst failure would initiate in the outer ring, whereas in case three it would most likely initiate in the 
middle ring, all else being equal. A similar trend is observed for the TMC cases, wherein case 5 (fig. 20(b)) would 
initiate burst failure in the outer ring and case 6 (fig. 20(c)) clearly would initiate failure in the inner ring given suffi- 
cient rotational speed. 

Since cases 2 and 5 had the most favorable material gradient with respect to minimizing the potential for disk 
separation (i.e., maximizing the speed of rotation prior to disk separation) the influence of misfit and rotational 
speed will be examined given these configurations. The results are presented in figures 21 to 23. In figures 21 and 
22 it is apparent that material gradation naturally leads to higher compressive (or at the least, smaller tensile) interfa- 
cial stress as the mean radius is increased until the ratio of thickness to mean radius becomes small and the influence 
saturates. Also, increasing the misfit (imposed hoop strain) leads as before to an increase in radial compressive 
stress. However, now since the outer disk is significantly stiffer than the inner disks, the compressive stress at inter- 
face A is now greater than I { . This result is in contrast to those given in figure 14. Further, for a fixed misfit strain, 
the radial stress generally increases with increasing mean radius; although this is dependent upon the magnitude of 
the mean radius. Given the previous discussion (see fig. 9) regarding temperature gradient results, one can envision 
that imposing an appropriate temperature gradient would lead to similar results as that of material gradation, since in 
essence the material properties as well as imposed internal thermal stresses will be nonuniform and naturally lead to 
increases in radial compressive stresses. 
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Finally, from figure 23, it can be clearly seen that increasing the rotational speed (CD), decreases the compressive 
stress at an interface (possibly eliminating it altogether depending upon the magnitude of mean radius, for a given 
misfit) as expected. One must keep in mind, that increasing the rotational speed significantly increases (proportional 
to or) the hoop (tangential) stress thus leading to potential burst failure of the system. Once again, from a design 
point of view all of these important parameters (misfit, mean radius, speed, material properties, and gradients of 
properties and/or loads) must be simultaneously optimized to achieve the “best” design. 


4.0 MATERIAL AND CONFIGURATION SELECTION 

In the previous section we spent a significant amount of time understanding the key design parameters that 
influence a single (or system of) rotating disks; such as mean radius, thickness, misfit, speed, material properties, 
and gradients of properties and/or loads. The interrelationship among these various parameters can be quite compli- 
cated and therefore typically requires a rigorous analysis be performed for each configuration (material and/or 
geometry) investigated. In this study we have focused primarily on two classes of composite materials, PMC and 
TMC, given a range (from relatively thick to thin) sizes of annular disks. In this section we will confine ourselves to 
single (solid or annular) disk configurations with constant height and attempt to determine the best configuration for 
maximizing energy storage. Clearly, an efficient flywheel stores as much energy per unit weight as possible, without 
“failing”. Failure can be defined in a number of different ways (e.g., uniaxial, multiaxial, burst criterion, or disk 
separation), induced as a result of several factors (e.g., over-speed (burst), decay of material properties (due to envi- 
ronmental factors), fatigue and/or creep conditions) and is dependent upon the specific application. Typically failure 
is dictated when the largest stress (i.e., hoop stress, see table VIII) due to centrifugal forces reaches the tensile 
strength (under burst conditions) or the fatigue strength (under cyclic conditions) of the material. The utilization, 
however, of anisotropic materials has made defining failure even more challenging as now the material strengths and 
properties are different in the various material directions and new modes of failure may be introduced. 

The “best” flywheel system is one that maximizes the kinetic energy per unit mass without failure. In practice, it 
should be remembered that the mass considered should be that of the whole system and the energy stored should be 
assessed only as the energy that can be supplied in normal service. Here, however, we will concern ourselves with 
the kinetic energy of the spinning disk, only. Consequently, given that this energy is defined as, U = Jc o 2 /2, with J 
being defined as the polar moment of inertia of the annular disk /Annular = n P f (b 4 - tf 4 )/2, or solid (with a = 0, i.e., 
/ so iid = rcp^ 4 /2) and the mass of the annular disk being w Annular = npt(b 2 - a 2 ), or solid (with a = 0, such that 
w solid = n P tb2 )i one can obtain easily the following ratio, which needs to be optimized, that is: 


or 



m 2 

(b 2 + a 2 \ 

= 0/ 


' m ' Annular 

4 J 


(38) 


(39) 


depending upon the desired geometry. 

To maximize the energy per unit mass, the task now becomes one of determining the maximum rotation speed 
(co) that a given disk configuration (be it material and geometric) can withstand before the induced stresses (or 
strains) exceed some allowable limit. In the case of determining the limit (burst) speed of the disk one would use the 
ultimate tensile strength of the material; whereas in the case of fatigue (where the disk is spun up and down 
repeatedly) one would most likely use the endurance limit as the design allowable stress. Note, that in the case of 
anisotropic materials (e.g., PMCs and TMCs) these allowable limits will be significantly different depending upon 
orientation and fiber volume fraction. For example in the case of PMCs the ratio of longitudinal to transverse ulti- 
mate tensile strength is approximately 30 to 1, whereas in the case of TMCs this ratio is 3 to 1. This together with 
the fact that the induced stress field is biaxial, makes the identification of an appropriate multiaxial failure criteria 
very important. 
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To illustrate this fact we will examine three separate failure conditions to determine the maximum speed of 
rotation for both a solid and annular flywheel. The first two conditions, limit the speed based solely on a single 
(uniaxial) stress component (i.e., tangential (longitudinal direction) and radial (transverse direction) stress, respec- 
tively). Alternatively, the third is based on an equivalent transversely isotropic effective (J 2 ) stress measure (multi- 
axial) described more fully in the next section, see equation (59), as well as in references 185 and 188. Returning to 
equations (36) and (37) it can be shown that the maximum rotation speed is either 
Hoop Stress Only 


or 

Radial Stress Only 



1 

b 2 f(a,b,$,r,v L ) 


or 

Multiaxial Condition 
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b 2 g(a,b,$,r,v L ) 
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where £, = — — 
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depending upon the desired failure condition utilized. The above expressions are quite general in that they can be 
used for both an annular ( a * 0) or solid ( a = 0) disk as well as a transversely isotropic (p < 1) or isotropic material 
(p =1). Ashby (ref. 182) has derived a similar expression to that of equation (40) for the special case of a solid, iso- 
tropic disk. In (ref. 182) Ashby denoted the term Ip as a performance index and employed it as a discriminator 
for the selection of the optimum material from which to construct a flywheel. The optimum material selected was a 
graphite (CFRP) or glass (GFRP) reinforced epoxy system, which can store between 150 to 350 KJ/kg of energy. It 
is interesting to note that the current work provides similar overall trends, however, the specific form now has an 
additional anisotropic factor that will modify the performance index suggested by Ashby (ref. 182) thereby making 
it more accurate. Inclusion of this additional factor requires that the appropriate anisotropic material properties be 
included in any database used to select between both monolithic (isotropic) and composite (anisotropic) materials. 

In the past others (e.g., see refs. (1) and (2)) have expressed the performance index coming from the specific 
energy for isotropic systems as 




m 


(43) 


where K was defined to be geometric shape factor. Furthermore, it was only noted in passing in the literature that 
this shape factor will become significantly more complex in the presence of anisotropic material behavior. Now by 
combining equations (39) and (42) one can define this complex “shape” factor, which accounts for both material 
anisotropy and stress state multiaxiality to be: 


K = 




y +/ 2 -,?-/} 


4b 2 


(44) 


Clearly this shape factor (only valid as shown here for constant height disks) is no longer just a geometric parameter 
but in actuality has both material and geometric attributes, K- K (P, alb, v L ) such that now both K and c /\ p 
become integral to the performance index. 
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Invoking the corresponding maximum speed resulting from the three criteria; the calculated stored kinetic 
energy per unit mass (specific energy), given an annular and solid disk made of the three PMC and TMC materials 
defined in tables VI and VII, are given in tables X and XI, respectively. Examining the results a number of conclu- 
sions become obvious. First, as expected from equations (38) and (40) to (42), the specific energy within a solid disk 
is independent of its geometry and only influenced by the degree of material anisotropy (stiffness and strength 
(UTS)). Secondly, for given annular disk, the specific energy increases as the mean radius is increased. This conclu- 
sion is particularly apparent when using the multiaxial criterion. Alternatively, however, the hoop only criterion 
does not show such a clear monotonic increase with increasing mean radius (as the location of the maximum hoop 
stress moves, this is particularly true for the case of strong anisotropy and relatively thick disks, see figs. 7 and 8), 
thus demonstrating the advantage of using a multiaxial criterion. Thirdly, it is apparent that the multiaxial criterion 
naturally accounts for the geometry of the disk; in that if a relatively thick anisotropic (with a weak bond between 
fiber and matrix) disk is analyzed the controlling stress is the radial component, whereas in the case of a relatively 
thin disk the hoop stress dominates. From tables X and XI it is apparent that the analyst would need to appreciate 
and account for this fact if they restricted themselves to using a simple uniaxial criterion. Fourthly, annular anisotro- 
pic disks are more efficient than solid anisotropic disks, see table XII. This is in sharp contrast to the fact that for 
isotropic materials, the specific energy of a solid disk is somewhere between 20 percent (for a thin annular disk) to 
100 percent (for a thick annular disk) more efficient than its annular counterpart. Note also, that if one where to only 
utilizes the hoop stress criterion, one would incorrectly conclude that in the case of only slightly anisotropic materi- 
als (TMCs) a solid disk is more efficient than an annular disk. Once again demonstrating the importance of using a 
multiaxial criterion in the presence of a multiaxial stress field. 

Finally, if we examine the ratio of the specific energy of PMC flywheels to that of TMC flywheels, it is appar- 
ent that flywheels made from PMCs are at least comparable to those made with TMCs, and more often they are 
approximately three times better (particularly for large sized flywheels (see table 13)). This result is not unexpected 
as PMCs are typically 2.5 times lighter than TMCs and as suggested by the above performance index, density plays 
a significant role. Also, as noted in table XIII the manufacturable fiber volume fraction of PMCs are typically sig- 
nificantly higher than TMCs and consequently the ultimate longitudinal strength of PMCs are typically significantly 
greater (1.5) than that of TMCs (compare tables XI and XII). However, the transverse strength of PMCs are 
approximately 0.2 that of TMCs; thus explaining why if one were to employ only a radial failure criterion, PMC 
flywheels would only perform approximately 0.3 to 0.5 as well as TMC flywheels, see table XIII. Remember how- 
ever, that thus far only a static failure condition (ultimate tensile strength) at room temperature has been used to 
determine the maximum allowable rotational speed; which is key to calculating the maximum stored energy. Under 
fatigue conditions and/or at elevated temperatures TMC flywheels may out perform flywheels made of PMCs. This 
will be a topic of future study. 


5.0 FAILURE ANALYSIS 

A key factor in any design is the ultimate load capacity of a given structure. Flywheels are no exception. The 
desire here is to maximize the energy storage per unit mass without experiencing catastrophic failure of the system. 
Again, for simplicity we will limit our discussion to that of a single rotating annular disk of constant height and 
examine the allowable limit speed under both monotonic and cyclic load histories. Multiple disk configurations with 
imposed misfit will be addressed in the future. 


5.1 Burst (Rupture) Limit 

As mentioned previously, numerous failure criterions can be defined, e.g., disk separation, delamination within 
a disk, or a burst criterion. Here our concern will be with calculating the burst pressure or burst speed of a single 
annular disk. Although, within the current RSL program disk separation is the primary design consideration, burst 
conditions must still be established. Currently, hydroburst tests are being considered as potential certification tests 
for establishing the quality of the flywheel configurations being manufactured. However, one should always remem- 
ber that for a general flywheel configuration the character of the problem changes significantly when going between 
internal pressure (hydroburst test) and rotation (spin test). That is, under internal pressure the ordering of the stresses 
is such that a r is always minimum (compressive), a. the intermediate and a 0 the maximum (tensile) value; whereas. 
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under rotation all stresses are tensile, with a r being either an intermediate value, in the case of plane stress, or a 
minimum in the case of plane strain. This stress state may be altered depending upon the design selected (e.g., a 
preloaded multidisk design will have compressive radial stress component throughout provided the preload is suffi- 
cient). Therefore, although spin tests are considerably more expensive than their hydroburst counterparts, they still 
remain the most prototypical and representative tests and thus the ultimate certification tests. Consequently, we are 
also interested in calculating the limit (burst) speed for a given rotating flywheel configuration. 

5.1.1 First Fracture . — An approach to calculating the burst limit for a circumferentially reinforced disk, for 
either internal pressure or rotation, is to define burst failure to be at the instant that the maximum hoop stress 
exceeds the material's ultimate tensile strength in the fiber direction (longitudinal). However, as we showed in the 
previous section, a more consistent criterion (than just the hoop stress alone) is to utilize an equivalent transversely 
isotropic effective stress (multiaxial) criterion. Consequently, the limit pressure and limit speed based on first frac- 
ture (at the location of maximum effective stress) are: 


and 


)limit 



1 


(45) 


® limit — 



(46) 


respectively. Where F.S. is the factor of safety imposed, g and/ were define previously, see equations (36) and (37) 
for the case of pure rotation, and now 



are taken directly from Arnold (ref. 175) for the case of internal pressure. The above first fracture based limits can 
be considered lower bounds on the corresponding limit loads, particularly for brittle like materials (e.g., PMC). 
However, in the case of more ductile materials (e.g., TMCs) one would prefer to employ plastic limit analysis to 
obtain a set of bounds. 
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5.1.2 Fullv-Plastic Limit . — Considering the case of the fully plastic, limit (burst) load (e.g., pressure or speed), 
extensive studies have been made previously for the case of isotropic tubes and disks under pressure and rotation; 
see Hill (ref. 183) and Nadai (ref. 184). Following these approaches, applications for the extended case of material 
transverse anisotropy under pressure for (plane-strain) cylinders and thin disks, were reported in Robinson and 
Pastor (ref. 185). A generalization of this anisotropic case under rotation will be outlined here. For convenience, we 
first summarize the governing equations of the present limit load problem, and outline the solution procedure. This 
will then be followed by a summary of the various proposals for the different limit loads. More detailed studies and 
derivations can be found in references 184 to 187. 

The required three conditions for the limit load solution are the satisfaction of (i) equilibrium (see eq. (1)) and 
any associated stress boundary conditions, (ii) compatibility for the “plastic" strain rates (see eq. (3)), and iii) the 
yield function (see eq. (59)). In particular when equilibrium and yield are combined we get for the generalized case 
of plane strain stress the following general equilibrium expression: 


where 



j ^4aK 2 +{b 2 -4ac}o? 
r 2 dr 
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with A, | and X 2 defined for the case of plane stress and plane strain as follows: 


Plane Stress: 
Plane Strain: 


— X«2 — 0 
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(50) 


Note the strength of anisotropy (Q is defined as, £= 4(^ 2 - l)/(4^ 2 - 1), where £ = a u /c u is the ratio of longitudi- 
nal to transverse yield (or ultimate) stress. Similarly, the required kinematic constrainf is tSe satisfaction of strain 
compatibility, that is: 


e r =-j-(r£e) (51) 

dr 

Finally, the yield (or ultimate) transverse shear stress, K , (first used in eq. (49)) is related to the longitudinal tensile 
yield stress by the following expression, 


2 _ (gj) 2 _ (g[) 2 
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(52) 
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Internal Pressure : As an example, for the case of internal pressure only, under plane stress, the final expression 
(including stress boundary conditions) for the lower bound of failure obtained from equation (49) will be: 


0 

i 

2K 


2dr\ 



(1 — T| 2 ) — t| 



^limit <- j 

2K 


(53) 


where generally this integration requires a numerical procedure. However, for the special case of plane strain under 
internal pressure only, analytical integration (between r = a, G r (a) = and r = b , G r (b) = 0) is possible leading 

to the following lower bound solution. 


flimit _ * fEf i n * 
2 K 2 V 1 a 


(54) 


It is also interesting to note the fact that the above solution also satisfies compatibility thus making it the “exact” 
limit pressure for the plane strain case and consequently a valid upper bound for the case of plane stress, as dis- 
cussed by Robinson and Pastor (ref. 185). 

Rotation : In the case of rotation (i.e., nonzero body-force term) the resulting differential equations for the lower 
bound estimates for both the cases of plane-strain and plane-stress are not separable and must be solved for numeri- 
cally. A straightforward iterative solution technique can be employed in conjunction with Euler integration until the 
boundary conditions, G r (a) = 0 and G r (b) = 0, are satisfied explicitly. Unfortunately, in the case of rotation the 
resulting stress state arising from the solution of equation (49) for the case of plane strain cannot be shown to be 
associated with a kinematically admissible strain field, as in the case of internal pressure loading. Consequently, 
both the plane stress and plane strain conditions merely yield lower bound (limit) solutions. 

An alternative analytical proposal for the upper, lower, as well as out-of-plane bounds have recently been 
developed and are described more fully in reference 187 with the final expressions taking on the following forms 
respectively: 

Lower Bound : 



(55) 


Upper Bound : 


Out of Plane Bound : 


where 
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Representative limit (upper and lower bound) curves for a given geometry (i.e., b/a = 1.8) are shown in 
figure 24; that is normalized effective pressure, PI 2K, versus strength of anisotropy, £. Note in the case of rotation 
this effective pressure is defined to be the “equivalent” centripetal load, P= pc o 2 (b 2 - a 2 )ll , and in the case of inter- 
nal pressure this is merely the applied internal pressure, P = P m . 

From figure 24 it is apparent that in the case of rotation the difference between upper and lower bounds is 
greater than in the case of internal pressure, for slightly anisotropic materials; with the averaged effective centripetal 
limit load (for an isotropic material, £= 0) being just slightly lower ( P/2K = 0.52) than the corresponding averaged 
limit pressure load ( P/2K « 0.58). Similarly, as the strength of anisotropy is increased (£> 0.8) it is interesting to 
note (and yet coincidental) that the plane strain solution for the case of internal pressure follows closely the average 
limit speed obtained for the rotation case. 

Also, from figure 24 we see that the critical value of reinforcement (referred to by Robinson and Pastor 
(ref. 185) as the condition of over-reinforcement) for the case of internal pressure occurs at a strength of anisotropy 
equal to, £ = 0.713 for a b/a = 1 .8. This point is denoted in figure 24 by the on-set of the plateau. Whereas in the 
case of rotation, a potential shear induced out-of-plane mode (see eq. (56)) is identified which is purely dependent 
upon the ratio of outer radius b , to height /?, (or length) of the disk (or cylinder). An insightful design oriented pre- 
sentation of this out-of-plane mode is shown in figure 25, where it is immediately apparent that the critical height (or 
length) to avoid (or cause) axial splitting before in-plane failure is dependent upon: (i) the strength of anisotropy, (ii) 
in-plane geometry of the disk {alb) and (iii) the ratio of in-plane to out of plane shear flow stress ( KIK *). 

These curves are simply obtained by equating a given in-plane bound to that of the out-of-plane bound. Though 
in actuality the modes will most likely interact in a complex way, following the notion of nonlinear interaction 
curves. Figure 25(a) corresponds to setting the lower in-plane bound equal to the out-of-plane bound while that of 
figure 25(b) corresponds to equating the average of the upper and lower in-plane bounds to that of the out-of-plane 
bound. Examining figure 25 it is clear that in the case of isotropy (£ = 0) the relative size of the disk (alb) has a sig- 
nificant impact on the allowable height (/?). Whereas, in the case of strongly anisotropic yielding (or failure) behav- 
ior, the relative size of the disk is not as important as is the ratio of the in-plane to out-of-plane shear flow stress (Kl 
K*). Furthermore, one can utilize the results displayed in figure 25 to either design a given disk to fail in-plane be- 
fore axial splitting takes place (i.e., adjusting h/b(K/K *) to be below the line) or to make sure that what is believed 
to be the more benign failure mode of axial splitting takes place prior to the more catastrophic in-plane failure (i.e., 
selecting an h/b(K/K*) such that it is above the associated line). 

Returning to the case of internal pressure, we see that the condition of over-reinforcement corresponds to a 
point in which the mode of failure switches from that of overall in-plane to that of localized out-of-plane flow within 
the vicinity of the inner radius (/* = a). This alternate upper bound limit is important as it suggests either the maxi- 
mum amount of reinforcement (volume fraction of fibers) one should utilize or alternatively the maximum relative 
size of the disk for a given amount of reinforcement. If one were to exceed either of these parameters, see figure 26, 
the resulting design would be inefficient and the resulting limit pressure would be limited by the value of the allow- 
able transverse shear flow stress of the material, as P = 2K. For example, in the case of PMC’s with relatively high 
volume fractions of fibers (£ > 0.95), the relative disk size, )JR m <0.15, should be small to ensure in-plane failure. 
From a practical point of view, it is interesting to note that from such a configuration one has the potential to obtain 
an accurate measure of the transverse shear flow stress ( K ), since at burst and for this geometry, P in = 2K. 

Recall, in section 4 we demonstrated that the stored energy per unit mass increases as the relative size ('klR fn ) of 
the disk is decreased. Consequently, one might conclude that the most efficient flywheel design is that of a relatively 
“thin” disk. If thicker disks are to be utilized, than one would naturally desire as high a transverse shear stress (K) as 
possible; thus suggesting the utilization of some type of 3-D weave so as to increase the transverse flow (ultimate) 
stress and consequently the controlling limit load of the given disk. 

Finally, it is important to realize that the first fracture analysis (using the effective J 2 criterion 1 ), initially dis- 
cussed, provides a minimum lower bound as compared to that coming from the ductile plastic limit analysis, see 
table XIV. This is consistent with the fact that in the first fracture case yielding (or fracture) is satisfied only at a 
single point, whereas in the plastic limit analysis yielding (or fracture, depending on the values used in the anisotro- 
pic strength measurement, Q is satisfied at all points throughout the disk. Consequently, one must be cognizant 
(when estimating the actual limit load) of whether the material tends to exhibit more ductile or brittle behavior dur- 


^ote that in the case of internal pressure a significant difference between the first fracture limit using maximum hoop stress only verses the 
effective J 2 criterion is observed (again since radial stress is always compressive the effective J 2 (multiaxial stress measure) is always larger 
than the hoop stress, maximum stress, alone). 
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in g failure, which is determined primarily by material selection (e.g. resign/sizing/flow) and the chosen manufactur- 
ing process (e.g., cure temperature, environment, etc.). In materials that exhibit both characteristics, interaction ef- 
fects between brittle and ductile failures become important, therefore this will be an area of future study. 


5.2 Fatigue Life 

Employing flywheels as energy storage systems inherently demand that the flywheel be repeatedly spun up and 
down, thus requiring the designer to be concerned about the fatigue capacity of any given flywheel system. To 
address this concern we will utilize a previously developed transversely isotropic fatigue damage model (refs. 188, 
189). This model has a single scalar damage state variable (damage is the same in all directions) that represents both 
mesoscale-initiation and -propagation of damage but whose rate of accumulation of damage is anisotropic (depen- 
dent upon the orientation of the preferred material direction relative to the loading direction). Furthermore, it has 
both a static fracture limit and endurance limit, accounts for mean stress effects and captures the typically observed 
nonlinear-cumulative effects in multiple block-program loading tests. The anisotropic fatigue (endurance) limit and 
static fracture surfaces are comprised of physically meaningful invariants that represent stress states that are likely 
to strongly influence the various damage modes within composites. For example, (i) the transverse shear stress 
(/j )— matrix cracking, (ii) the longitudinal shear stress (A,)- dictates interfacial degradation, and (iii) the maximum 
normal stress (/ 3 ) in the fiber direction which dictates fiber breakage. All three of these invariants are combined into 
a function representing the effective transversely isotropic J 2 invariant: 




!(), 


(4£ n 2 -l) 9 

(4^ ) — 1 )/ 1 H =; /-> H /-< 
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(59) 


where 

h=J 2 -/+L 3 



in which 


J 2 S ij S ij 
I = DjjSjj 
I = DjjSj k S kj 
^ij — djdj 

$ij = ®ij ~ ^ ®kfiij 

and d; (/ = 1,2,3) are the components of a unit vector denoting the local fiber direction. Note this effective J 2 invari- 
ant (eq. (59)) has been the multiaxial criterion used throughout this study, but specialized to the pertinent case of 
circumferential reinforcement with rj = 1 (equal longitudinal and transverse shear response). This special case is 
simply written in terms of the radial and hoop stress as follows: 


F( } = )ct " +a ® _a,ae } (60) 

where the open bracket subscript notation in equation (59) and the above is used to distinguish between the static 
fracture, endurance limit and normalization surfaces, see reference 189. 
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This fatigue model has already been successfully applied to the problem of a circumferentially reinforced disk 
subjected to internal pressure, see (refs. 190 and 191 ). Consequently, here we will focus our attention only on the 
fatigue life of a circumferentially reinforced, rotating disk and in particular on the influence of the key geometric 
parameters, i.e., mean radius and thickness. The fatigue life model of (ref. 188) requires at a minimum the availabil- 
ity of both longitudinal and transverse S-N curves to describe the evolution of damage for a given material system. 
As this data is presently not available for the PMC system of interest we will confine our investigation to a previ- 
ously examined SiC/Ti system; the associated damage model parameters being completely defined in reference 
(ref. 191). 

Results for the case of a single, matrix cladded (inner and outer monolithic disk), disk subjected to multiple 
cycles of rotational loading (triangular waveform) is shown in figures 27 to 29. In figure 27 the limit burst speed 
versus the geometry of the disk is presented, wherein this limit speed corresponds to a given induced stress state at 
the various radial locations (i.e., inner radius (ID) and outer radius (OD)) that fulfills the static fracture condition of 
equation (51), that is F L - 1 = 0. This analysis is defined as an uncoupled approach, see reference 190. Figure 27(a) 
clearly shows that the limit speed greatly decreases with increasing mean radius; while figure 27(b) illustrates that as 
the relative thickness (k!R m ) of the disk compared to the mean radius increases the difference between the inner and 
outer limit speed increases greatly as well, as one would expect. In figure 28, the resulting speed versus cycles to 
failure curve defining fatigue initiation at the inner radius (ID) of the composite core for various mean radii are dis- 
played. Again, we observe the expected trend associated with increasing the mean radius of the disk; that is, a 
marked decrease in either the limit speed (for similar life times) or number of cycles to end of life (for the same 
speed). If one where to normalize the calculated fatigue limit speed with the corresponding static limit (burst) speed 
at the same radial location for the given geometry one can obtain a master fatigue initiation life curve that is essen- 
tially independent of mean radius variation, see figure 29. A similar master curve can be obtained corresponding to 
fatigue initiation at other radial locations within the disk as well. Also, if one were to conduct a coupled deformation 
and fatigue damage analysis, as done previously for the case of internal pressure (see refs. 190 and 191) a single 
master design life curve representing the global failure of the disk due to fatigue could be produced. Such a curve 
will be produced in the near future for the PMC reinforced disks of interest to the RSL program when the necessary 
experimental data for characterizing the current life model become available. 


6.0 CONCLUSIONS 

In this study we have attempted to put into perspective the problem of a rotating disk, be it a single disk or a 
number of concentric disks forming a unit. An analytical model capable of performing an elastic stress analysis for 
single/multiple, annular/solid, anisotropic/isotropic disk systems, subjected to both pressure surface tractions, body 
forces (in the form of temperature-changes and rotation fields) and interfacial misfits has been derived and dis- 
cussed. Results of an extensive parametric study were presented to clearly define the key design variables and their 
associated influence. In general the important parameters are misfit, mean radius, thickness, material property and/or 
load gradation, and speed; all of which must be simultaneously optimized to achieve the “best” design. All of the 
observations made herein regarding the deformation analysis are in total agreement with previous conclusions and 
observations made in the past by other investigators. 

Also, the important issue of defining proper performance/merit indices (based on the specific stored energy), in 
the presence of multiaxiality and material anisotropy have been addressed and utilized to discuss the difference 
between flywheels made from PMC and TMC materials with either an annular or solid geometry. An interesting 
observation made is that an annular anisotropic disk is more efficient than its solid counterpart; this is in sharp con- 
trast to the fact that an isotropic solid disk is more efficient than an annular disk. 

Finally two major aspects of failure analysis, that is the static and cyclic limit (burst) speeds were addressed. In 
the case of static limit loads, upper, lower, and out-of-plane bounds for disks with constant thickness were presented 
for both the case of internal pressure loading (as one would see in a hydroburst test) and pure rotation (as in the case 
of a free spinning disk). The results (interaction diagrams) were displayed in a designer friendly format. For the case 
of fatigue, a representative fatigue/life master curve was illustrated in which the normalized limit speed versus num- 
ber of applied cycles was described for a cladded TMC disk application. 
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TABLE I.— ENERGY STORAGE TYPES (ref. 5) 


Storage type 

Specific energy, 
Wh/kg 

Magnetic fields 


Superconducting coil 

1-2 

Elastic deformations: 


Steel spring 

0.09 

Natural rubber band 

8.8 

Electrochemical reaction 


Lead-acid battery 

17.9 

Nickel-cadmium battery 

30.6 

Kinetic energy 


Maraging steel flywheel 

55.5 

4340 steel flywheel 

33.3 

Composite flywheel 

a 213.8 


a Using longitudinal strength of Kevlar 


without any shape factor. 


TABLE II.— ROTATING DISKS COMPSED OF ISOTROPIC MATERIALS 


Conditions/type of study 

Analytical 

Numerical 

Experimental 

Variable thickness 

[5]. [6], [7], [16], [17], [20], 
[22], [23], [26], [32], [44], 
[55], 

[8], [12], [24], [30], 
[42], [49], [50], 

[28], [31], [40], 

Uniform thickness 

[5], [6], [7], [16], [17]. [18], 
[19], [20], [21], [22], [23], 
[25], [26], [32], [33], [36], 
[44], [45], [51], [55], 

[8], [12], [24], [27], 
[30], [34], [35], [37], 
[38], [39], [41], [42], 
[43], [46], [47], [49], 
[50], [51], [53], [54], 

[28], [31], [40], 

Misfit 


m. 


Temp 

[18], [48], 

[39], [46], 

[31], 


TABLE III.— ROTATING DISKS COMPOSED OF ANISOTROPIC MATERIALS 


Conditions/type of study 

Analytical 

Numerical 

Experimental 

Variable thickness 

[59], [61], [64], [66], [71], [88], 

[9], [10], [107], [109], [116], 


Uniform thickness 

[5], [56], [57], [58], [59], [60], 
[61], [62],[63], [64], [66], [67], 
[68], [69], [70], [71], [73], [75], 
[76], [88], [92], [95], [101], 
[103], [113], 

[65], [72], [74], [77], [82], 
[84], [94], [98], [100], [104], 
[105]. [106], [107], [108], 
[109], [1 10], [1 14], [116], 

[60], [74], [75], [78], [79], 
[80], [81], [83], [85], [86], 
[87], [89], [90], [91], [93], 
[96], [97], [99], [102], 
[110], [112], [114], 

Misfit 

[5], [70], 

[94], [111], 

[90], 

Temp 

[61], [67], 

[74], [82], [100], [108], 

[74], [99]. 
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TABLE IV.— NON-LINEAR ANALYSIS CLASSIFICATION FOR ROTATING DISKS 



Conditions/type of study 
Elasto-Plastic 
Variable thickness 


Isotropic 


Anisotropic 


Creep/Viscoelastic 
Variable thickness 


TABLE V.— SUMMARY OF PROBLEMS INVESTIGATED TO VALIDATE 
PRESENT ANALYTICAL SOLUTION. 


Type/Altemate Method 

Rotation 

Pressure 

The rmal 

Misfit 

Solid Anisotropic Disk 
FEA 

X 

X 

X 


Single Annular Anisotropic Disk 
Lekhnitskii [58] & FEA 

X 




Arnold [175] 


X 

X 


Hub/ Anisotropic Disk 
FEA 



X 

X 

Two Annular, isotropic Disks 
Ugural and Fenster [176] 




X 


TABLE VI.— ROOM TEMPERATURE PMC MATERIAL PROPERTY DATA 


Property 

Fiber, 
(ref. 180) 

Resin, 
(ref. 181) 

v f = 40 percent 

v f = 60 percent 

v f = 80 percent 

£,, Msi 

40.0 

0.677 

16.4 

23.1 

32.1 

£ y, Msi 

2.0 

0.677 

1.06 

1.26 

1.58 

V L 

0.25 

0.33 

0.295 

0.282 

0.267 

v r 

0.25 

0.33 

0.42 

0.373 

0.303 

G Msi 

2.9 

0.2544 

0.425 

0.62 

1.5 

a L , in./in./°F 

-0.22XKT 06 

42.0X10- 06 

0.845X10" 06 

0.269X10- 06 

-O.O32xlO r06 

a r , in./in./°F 

6.28X10- 06 

42.0x1 O' 06 

33.7x10^ 

25.4x10^ 

17.0XKT 06 

p, kip-sec 2 /in 4 

1.651 lxlO- 07 

1.2163XKT 07 

1.3033XKT 07 

1.477xl<r° 7 

1.5657XKT 07 

UTS, Ksi 

747 to 790 

17.5 

a l 10 

a 302 L 
b 10T 

a 398 


a Determined via rule of mixture with a 0.3 knockdown for in-situ strength. 
Suggested from experiment. 
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TABLE VII.— ROOM TEMPERATURE TMC MATERIAL PROPERTY DATA 


Property 

Fiber 
SCS-6, 
(ref. 192) 

Matrix 
Ti-15-3, 
(ref. 192) 

SCS6/Ti-15-3 
v f =15 percent 

SCS6/Ti-15-3 
v f =35 percent, 
(ref. 192) 

SCS6/Ti-15-3 
v f =55% 

E,, Msi 

58 

13.0 

19.8 

28.8 

26.1 

37.8 

E 7 . Msi 

58 

13.0 

15.7 

20.2 

17.4 

27.0 

v r 

0.19 

0.33 

0.306 

0.272 


0.249 

v r 

0.19 

0.33 

0.366 

0.37 


0.319 

G/, Msi 

24.36 

4.924 

5.65 

7.0 


9.62 

a L , in/in/°C 

2.2x1 (P 06 

8.IxlO _0t> 

5.54x10-°'’ 

3.99X10- 06 


3.18 xlO- 06 

a T , in/in/°C 

2.2XKT 06 

8-lxlCT 06 

7.66x10^ 

6.38x10-°° 


5.24 xlO- 06 

p, kip-sec 2 /in 4 

3.085X10- 07 

4.21xl0^ 7 

4.04x1 0 -07 

3.81x10-°’ 


3.591X10- 07 

UTS, Ksi 

550 to 650 

126.2 

a 185 

a 257 

201.6 L 
60.9 T 

a 347 


“Determined via rule of mixture: Bolded values are experimentally measured, (ref. 192). 


TABLE VIII.— MAXIMUM NORMALIZED RADIAL AND TANGENTIAL STRESS FIELDS 


R m 

o r /pco : =g( (3. R m . X) b ; 

<ypar=f(P. R m . X) b : 

X=2 

'Ti- 

ll 

it 

00 

II 

K) 

n 

00 

II 


PMC 

TMC 

PMC 

TMC 

PMC 

TMC 

PMC 

TMC 

PMC 

TMC 

PMC 

TMC 

4 

1.08 

1.6 

2.1 

6.3 

3.8 

19.5 

17.8 

22.0 

22. 

29. 

39. 

25 

8 

1.46 

1.64 

4.3 

6.5 

8.5 

25.1 

72.8 

75.3 

71 

88 

88 

116 

16 

1.6 

1.65 

5.8 

6.6 

17.4 

26.1 

276 

278 

291 

301 

286 

353 

32 

1.64 

1.65 

6.4 

6.5 

23.4 

26.3 

1066 

1067 

1106 

1111 

1165 

1205 


TABLE IX.— DESCRIPTION OF MATERIAL GRADATION BY CASE 


Case 

Disk 1 

Disk 2 

Disk 3 

One 

v f = 60%, PMC 

v f = 60%, PMC 

v f = 60%, PMC 

Two 

v,= 20%. PMC 

v f = 60%. PMC 

v f = 80%, PMC 

Three 

v t = 80%. PMC 

v f = 60%. PMC 

Vf = 20%. PMC 

Four 

v f = 35%, SCS-6/Ti-51-3 

v f = 35%, SCS-6/Ti-51-3 

Vf = 35%, SCS-6/Ti-51-3 

Five 

v f = 15%, SCS-6/Ti-51-3 

v f = 35%. SCS-6/Ti-51-3 

Vf = 55%, SCS-6/Ti-51-3 

Six 

v f = 55%, SCS-6/Ti-51-3 

v,= 35%. SCS-6/Ti-51-3 

v f = 15%, SCS-6/Ti-51-3 


TABLE X.— MAXIMUM SPECIFIC STORED ENERGY (KJ/KG) FOR A SINGLE 
DISK MADE WITH PMC 


Criterion 

Annular Disk Case ( X = 3) 

Solid Disk Case (a = 0) 


^=4“ | R m = 8 | R m = 16 | R m = 32 

b = 5.5 | b = 9.5 | b = 17.5 | b = 33.5 

Hoop Only 

Vf = 40% 

327.86 

295.29 

299.96 

313.04 

275.62 

275.67 

275.74 

275.79 

Vf = 60% 

792.16 

713.46 

724.75 

756.37 

665.94 

666.06 

666.23 

666.36 

Vf = 80% 

983.93 

886.17 

900.19 

939.46 

827.14 

827.29 

827.51 

827.67 

Radial Only 

Vf = 40% 

113.74 

250.12 

825.3 

3127.1 

87.285 

87.357 

87.514 

87.64 

Vf = 60% 

267.56 

588.40 

1942.1 

7357.1 

205.34 

205.51 

205.88 

206.17 

v f = 80% 

338.8 

745.05 

2458.1 

9316.1 

260.01 

260.22 

260.69 

261.06 

Multiaxial 

Vf = 40% 

115.36 

211.14 

299.96 

313.04 

89.848 

90.048 

90.109 

90.01 

Vf = 60% 

278.72 

510.15 

724.75 

756.37 

217.09 

217.57 

217.72 

217.48 

Vf = 80% 

346.19 

633.65 

900.19 

939.46 

269.64 

270.24 

270.42 

270.12 


a R m , b and X has the units of inches. 
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TABLE XI.— MAXIMUM SPECIFIC STORED ENERGY (KJ/KG) FOR A SINGLE 
DISK MADE WITH TMC 


Criterion 

Annular Disk Case ( X = 3) 

Solid Disk Case (a = 0) 


R m = 4°| = 8 1 R„,= 16 I R„T32~ 

b = 5.5 | b = 9.5 | b = 17.5 | b = 33.5 

Hoop Only 

v f = 15% 

105.53 

119.76 

131.69 

139.10 

190.68 

190.75 

190.68 

190.69 

v f = 35% 

155.24 

176.18 

193.74 

204.63 

280.51 

280.61 

280.50 

280.53 

v,= 55% 

222.74 

252.79 

277.97 

293.60 

402.48 

402.62 

402.46 

402.50 

Radial Only 

v f = 15% 

225.36 

802.10 

3119.1 

12361.0 

72.94 

72.97 

72.94 

72.97 

v f = 35% 

332.41 

1183.10 

4610.0 

18231.0 

107.58 

107.64 

107.58 

107.63 

v f = 55% 

452.87 

1612.10 

6267.1 

24831.0 

146.56 

146.64 

146.57 

146.63 

Multiaxial 

Vf = 15% 

105.53 

1 19.76 

131.69 

139.1 

71.79 

71.8 

71.79 

71.84 

Vf := 35% 

155.24 

176.18 

193.74 

204.63 

105.61 

105.62 

105.62 

105.69 

Vf= 55% 

222.74 

252.79 

277.97 

293.6 

151.53 

151.54 

151.54 

151.64 


a R m , b and X has the units of inches. 


TABLE XII.— RATIO OF ANNULAR DISK TO SOLID DISK MAXIMUM SPECIFIED STORED 

ENERGY 


m hu\ 

PMC (X = 3) 

TMC (X = 3) 

\m) A l 

R m = 4 a 

R m = 8 

R m =16 

R m = 32 

R m= 4 

PO 

3 

II 

00 

R m = 16 

<N 

CO 

II 

E 

PC 

Criterion 

Hoop Only 

1.19 

1.071 

1.088 

1.135 

0.553 

0.628 

0.691 

0.729 

Radial Only 

1.303 

2.863 

9.43 

35.686 

3.09 

10.991 

42.758 

169.359 

Multiaxial 

1.284 

2.345 

3.329 

3.478 

1.47 

1.668 

1.834 

1.936 


a R m and Xhas the units of inches. 


TABLE XIII.— RATIO PMC TO THAT OF TMC 
FLYWHEELS SPECIFIC ENERGY 


(-) /n 

V m ) PMC ! V m )tMC 

Annular Disk Cases (X = 3) 

- 3 - 

n 

E 

PC 

po 

3 

II 

00 

R m =16 

Rm = 32 

Hoop Only 

40%/ 15% 

3.107 

2.466 

2.278 

2.25 

60%/35% 

5.103 

4.05 

3.741 

3.696 

80%/55% 

4.417 

3.506 

3.238 

3.2 

Radial Only 

40%/ 15% 

0.505 

0.312 

0.265 

0.253 

60%/35% 

0.805 

0.497 

0.422 

0.404 

80%/55% 

0.748 

0.462 

0.392 

0.375 

Multiaxial 

40%/15% 

1.093 

1.763 

2.278 

2.25 

60%/35% 

1.795 

2.896 

3.741 

3.696 

80%/55% 

1.554 

2.507 

3.238 

3.2 


a R m and Xhas the units of inches. 


TABLE XIV.— COMPARISON OF BRITTLE AND 
DUCTILE LIMITS 


b/a = 1.35, £=0.93 

Limit pressure, 
Ksi 

Limit speed, 
rpm 

Ductile 

— 

— 

Upper Bound 

61.038 

55 320 

Lower Bound 

54.237 

54 760 

Brittle 



First Fracture (/, ) 

39.1 

51 180 

First Fracture (o max ) 

58.3 

51 180 
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Fibers 
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Figure 1 . — Single disk, reinforced in the circumferential 
direction. 



L, e 


22, T, r 



Figure 2. — Schematic of "n" concentric disks. 
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Figure 3. — Comparison between the current study and 
Arnold [175]. 




Figure 4. — Comparison between the current study and 
Lekhnitskii [58]. 
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0.2 



r 
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Figure 5. — Normalized radial and tangential stresses. 
Comparison with maximum given by Ugural and 
Fenster [1 76]. 



Misfit 


Figure 6. — Schematic representation of two concentric 
disks with misfit. 
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PMC , X = 3 


PMC , X = 3 



(c) R m W R m 

Figure 7. — Normalized radial “g” and tangential “f” stresses (aka distribution factors) versus normalized radius 
given a single rotating disk made from PMC. Figure (a) and (b) correspond to a 3 in. thick disk, while (c) and (d) 
to a 7 in. thick disk. 
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(a) 


(b) 


TMC, X = 7 



(c) R m 



(d) R m 


Figure 8. — Normalized radial “g” and tangential “f” stresses (aka distribution factors) versus normalized radius 
given a single rotating disk made from TMC. Figure (a) and (b) correspond to a 3 in. thick disk, while (c) and (d) 
to a 7 in. thick disk. 


NASA/TM — 200 1 -210578 


40 



Applied 

temperature 

profile 


T-T 0 

profile 


'—Temperature 


Constant 

Constant 

Decrease 

1 

Increase 


i t 

) a t 

■ Constant temperature 


— — ■ Linearly increasing from 
inside temperature 




r 


Constant temperature 


— ■ Linearly increasing from 
inside temperature 
Linearly decreasing from 
inside temperature 



-20 


0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 


r 



A, 

in. 



Figure 9. — Influence of temperature gradient on a 
single disk. 


Figure 10. — Influence of applying pressure along the 
inner surface of a single disk. A - thickness of disk. 
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Uniform temperature 
— — Linearly increasing from 
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■ ■ Uniform temperature 
— — Linearly increasing from 



r 


2.0 - 
1.5 — 


■ - Uniform temperature 

“ — Linearly increasing from 
inside temperature 
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0.0 1 1 1 1 1 1 1 1 

0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 
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Figure 11. — Influence of thermal, internal pressure, 
and rotation on a single disk. Speed of rotation is 
15 000 rpm. 


Figure 12. — Influence of thermal, internal pressure, 
and rotation on a single disk. Speed of rotation is 
60 000 rpm. 
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Misfit 8 2 — \ 

\ 

\ 


Figure 1 3. — Schematic representation of three con- 
centric disks with misfit. 



(C) R m 

Figure 14. — Normalized radial stress for disk sub- 
jected to misfit and rotation. Case (a) has no misfit, 
i. e., e h = 0, (b) has £ h = 0.1 % and (c) has E h = 0.2% 
applied at all interfaces. 
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3.2 i— 



Figure 15. — Normalized tangential stress for disk sub- 
jected to misfit and rotation. Case (a) has no misfit, 
i. e., e h = 0, (b) has e h = 0.1 % and (c) has e h = 0.2% 
applied at all interfaces. 


Figure 16. — Influence of varying internal disks 
thickness on the radial and tangential stresses. 

{2, 0.5, 0.5}, X 2 {0.5, 0.5, 2}, X 3 {0.5, 2, 0.5}. The 
misfit here is £ h = 0.1 %. 
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2 



in. 




Figure 17. — Radial stress for different material 
configurations subjected to rotational speed 
of 30 000 rpm. (a) material case 1 , (b) ma- 
terial case 2 and (c) material case 3. (see 
Table 9) 



(a) r, in 




Figure 18.— Tangential stress for different ma- 
terial configurations subjected to rotational 
speed of 30 000 rpm. (a) material case 1 , 

(b) material case 2 and (c) material case 3. 
(see Table 9). 
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200 



Figure 19. — Radial stress for different material 
configurations subjected to rotational speed 
of 30 000 rpm. (a) material case 4, (b) material 
case 5 and (c) material case 6. (see Table 9). 


Figure 20. — Tangential stress for different material 
configurations subjected to rotational speed of 
30 000 rpm. (a) material case 4, (b) material 
case 5 and (c) material case 6. (see Table 9). 
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Figure 21 . — Radial stress for material case 2 for 
different applied misfits subjected to rotational 
speed of 30 000 rpm. (a) e h = 0, (b) e h = 0.1 % 
and (c) £ h = 0.2%. 


Figure 22. — Radial stress for material case 5 for 
different applied misfits subjected to rotational 
speed of 30 000 rpm. (a) £ h = 0, (b) £ h = 0.1 % 
and (c) £ h = 0.2%. 
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3.0 



(a) r, in 



(b) r, in 



(c) r, in 


Figure 23. — Radial stress for material case 2 for 
different speeds of rotation, (a) to = 15 000 rpm, 
(b) co = 30 000 rpm, and (c) co = 60 000 rpm. 



Strength of anisotropy, £ 

Figure 24. — Illustrates the upper and lower bounds 
for both rotation (solid lines) and internal pressure 
(symbols) limit loads by plotting the normalized 
effective pressure versus strength of anisotropy. 




0.0 0.2 0.4 0.6 0.8 1.0 

(b) £ 

Figure 25.— Interaction curves for the out-of-plane 
failure bound: (a) represents the limit as com- 
pared to the in-plane lower bound, and (b) 
represents the limit as compared to the average 
in-plane bound. The h is the height of disk, b outer 
radius, £ the strength of anisotropy. 
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£ 


Figure 26. — Illustrates the interaction curve differ- 
entiating in-plane failure (below the line) and 
out-of-plane flow (on or above the line). | is the 
ratio of outer to inner radius and £ is a measure of 
the strength of anisotropy. 


E 

Q. 

h 

3 



(a) Mean radius affect 


E 

Q. 

ll 

3 



(b) Thickness affect 


Figure 27. — Static limit or burst speed as a function 
of geometry, (a) maximum speed versus mean 
radius, (b) maximum speed versus the ratio of disk 
thickness to that of mean radius. 
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Fiber 



Figure 28. — Generalized speed versus cycles curve indicating the onset 
of fatigue initiation cracks at the inner radius of the composite core. 



Fatigue life, N f 

Figure 29. — Normalized speed versus cycle master curve denoting 
fatigue initiation at the inner diameter. to L is the static limit speed 
and o> is the maximum applied speed within a cycle. 
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